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Using 448.0 x 10° w 3686) events collected with the BESIII detector, an amplitude analysis is 
performed for (3686) > yxe1, Xeı  nata decays. The most dominant two-body structure 


observed is ao(980)*2*; ao(980)* — nr. 


The ao(980) line shape is modeled using a dispersion 


relation, and a significant nonzero ao(980) coupling to the n'm channel is measured. We observe 
Xe1 — a2(1700)7 production for the first time, with a significance larger than 170. The production 
of mesons with exotic quantum numbers, JPO = 17+, is investigated, and upper limits for the 


branching fractions x21 — 71(1400) 


PACS numbers: 13.25.Gv, 14.40.Be, 14.40.Pq, 14.40.Rt 


I. INTRODUCTION 


Charmonium decays provide a rich laboratory for light 
meson spectroscopy. Large samples of charmonium states 
with JPO = 177, like the J/w and (3686), are easily 
produced at e*e7 colliders, and their transitions pro- 
vide sizable charmonium samples with other JPO quan- 
tum numbers, like the xa (1717). The yea — Tum 
decay is suitable for studying the production of exotic 
mesons with JPO = 17+, which could be observed de- 
caying into the mz final state. The lowest orbital ex- 
citation of a two-body combination in xa decays to 
three pseudoscalars, for instance xa — nmr, is the S- 
wave transition, in which if a resonance is produced, 
it has to have JPO = 1-+. Several candidates with 
JPC = 17+, decaying into different final states, such as 
nn, nm, £1(1270)7, bi(1235)z and pr, have been reported 
by various experiments, and these have been thoroughly 
reviewed in Ref. [1]. The lightest exotic meson candidate 
is the 7,(1400) |2], reported only in the nz final state 
by GAMS [3], KEK (4), Crystal Barrel [5], and E852 [6], 
but its resonance nature is controversial [7]. The most 
promising JPO = 17+ candidate, the zi (1600) [2], could 
also couple to the 77, since it has been observed in the 
n'r channel by VES [8] and E852 [g). 


EnF, Xe — mı(1600)* rF, and xe — mı(2015)*rF, with 
subsequent mı ( X)” — mz decay, are determined. 


The CLEO-c collaboration reported evidence of an ex- 
otic signal in ye. — n'rtn—_ decays, consistent with 
(1600) — n'm production [10]. However, other pos- 
sible exotic signals that could be expected have not been 
observed in either Xea 3 nrn or Xea — at de- 
cays. With a more than 15 times larger data sample at 
BESIII, there is an opportunity to search for the pro- 
duction of mı exotic mesons. In this work we investi- 
gate possible production of exotic mesons in the mass 
region (1.3-2.0) GeV/c’, decaying into the mrt + c.c. fi- 
nal state, namely the 7(1400), 71(1600), and (2015), 
using Yer > natn decays. Charge conjugation and 
isospin symmetry are assumed in this analysis. 


Additional motivation for studying these decays is that 
a very prominent a9(980) — nz signal of high purity 
was observed in xa — nwt, by CLEO-c . The 
ao(980) was discovered several decades ago, but its na- 
ture was puzzling from the beginning, leading to the 
hypothesis that it is a four-quark rather than an ordi- 
nary qq state aHa. The first coupled meson-meson 
(mm, KK, n'r) scattering amplitudes based on lattice 
QCD calculations indicate that the ao(980) might 
be a resonance strongly coupled to nr and K K channels, 
which does not manifest itself as a symmetric bump in 
the spectra. Recent theoretical work based on the chiral 
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unitarity approach also points that the ag(980), as well 
as the ø and fo(980) states, could be dynamically gener- 
ated through meson-meson interactions, for example in 
heavy-meson decays: Xe1 — NTT and Ne > NTT [16]. 
However, there is still no consensus on the exact role that 
meson-meson loops play in forming of the ao (980), which 
is now generally accepted as a four-quark object, see 
and reference therein. 


The ao(980) indeed decays dominantly into nz and 
KK final states; the latter has a profound influence on 
the ag (980) line shape in the 77 channel, due to the prox- 
imity of the KK threshold to the ao(980) mass. Differ- 
ent experiments, E852 KEJ Crystal Barrel ig, and 
CLEO-c analyzed data to determine the couplings of 
the ao(980) to the nT (gnr) and KK final states (ggg), 
in order to help resolve the true nature of the ao(980). 
This is not an exhaustive list of analyses: it points out 
that the values obtained for the ao (980) parameters vary 
considerably among various analyses. 


Another channel of interest is a9(980) — wm, with 
the threshold more than 100 MeV/c? above the ag(980) 
mass. The first direct observation of the decay ao (980) > 
n'r was reported by CLEO-c Ed, using a sample of 
26 x 10° (3686) decays. The ao(980) coupling to the 
n'm channel, gyn, was determined from Xe — nnn 
decays, although the analysis was not very sensitive to 
the ao(980) — n'r component in the a9(980) — nz in- 
variant mass distribution, and g,, was found to be con- 
sistent with zero. In many analyses of ag(980) couplings, 
gn'n has not been measured. For example, its value was 
fixed in Ref. based on SU(8) flavor-mixing predic- 
tions. Using a clean sample of xai produced in the radia- 
tive transition w(3686) > yxa at BESIII, we investigate 
the xa > natn decays to test if the a9(980) — n7 in- 
variant mass distribution is sensitive to 7' production. 
Dispersion integrals in the description of the a9(980) line 
shape are used to determine the ao(980) parameters, its 
invariant mass, ™,,(9g0), and three coupling constants, 
Jnn, IKK and gyr. This information might help in de- 
termining the quark structure of the ao(980). 


In this xai decay mode, it is also possible to study 
Xe. — @2(1700)7; a2(1700) — nr production. The 
a2(1700) has been reported in this decay mode by Crys- 
tal Barrel and Belle Rd, but still is not accepted 
as an established resonance by the Particle Data Group 
(PDG) [2]. 


II. EVENT SELECTION 


events, collected in 2009 and 2012 with the BE- 
SIII detector (25). We select 95% of possible 7 decays, in 
the n — yy, n > at z and n — 7°r°r® decay modes. 
For each 7(3686) > ynrtn final state topology, exclu- 
sive Monte Carlo (MC) samples are generated according 


For our studies we = = Sor 10° (3686) 


to the relative branching fractions given in Table[I] equiv- 
alent to a total of 2 x 107 7)(3686) > YXc1; Xe — NT ta 
events. The background is studied using an inclusive MC 
sample of 106x 10° generic (3686) events. 


BESIII is a conventional solenoidal magnet detector 
that has almost full geometrical acceptance, and four 
main components: the main drift chamber (MDC), elec- 
tromagnetic calorimeter (EMC), time-of-flight detector, 
all enclosed in 1 T magnetic field, and the muon cham- 
ber. The momentum resolution for majority of charged 
particles is better than 0.5%. The energy resolution for 
1.0 GeV photons in the barrel (end-cap) region of the 
EMC is 2.5% (5%). For the majority of photons in the 
barrel region, with the energy between 100 and 200 MeV, 
the energy resolution is better than 4%. Details of the 
BESIII detector and its performance can be found in 
Ref (25). 


Good photon candidates are selected from isolated 
EMC showers with energy larger than 25 (50) MeV in 
the barrel (end-cap) region, corresponding to the polar 
angle, 0, satisfying | cos 0| < 0.80 (0.86 < | cos 0| < 0.92). 
The timing of good EMC showers is required to be within 
700 ns of the trigger time. Charged tracks must satisfy 
|cosð| < 0.93, and the point of closest approach of a 
track from the interaction point along the beam direc- 
tion is required to be within 20 cm and within 2 cm 
perpendicular to the beam direction. All charged tracks 
are assumed to be pions, and the inclusive MC sample is 
used to verify that the kaon contamination in the final 
sample is negligible in each of the 7 channels. We require 
two charged tracks for the 7 — yy and n — 37° chan- 
nels, and four tracks for the n + ntr ar? channel, with 
zero net charge. For n — yy and n —> m+2~7°, at least 
three photon candidates are required, and for n — 37° 
at least seven photon candidates. The invariant mass of 
two-photon combinations is kinematically constrained to 
the zŠ or 7 mass. 


The sum of momenta of all final-state particles, for 
a given final state topology, is constrained to the ini- 
tial 7)(3686) momentum. If multiple combinations for an 
event are found, the one with the smallest NGG is re- 
tained. Here NC refers to the number of constraints, 
which is four plus the number of two-photon 7° and n 
candidates in the final state (see Table[I). 


A. Xa > nrin event selection 


The xa. > mm! candidates in 7 three-pion decays 
are selected by requiring that the invariant mass of three 
pions satisfy 


0.535 < m(3m) < 0.560 GeV /2. (1) 


For the ņn — yy candidates, we require that the mass 
constraint fit for n — yy satisfies x2, < 15. The 
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TABLE I. Characteristics of the ņ decay channels used to 
reconstruct the (3686) — ynatm~ decays: branching frac- 
tion B, final state topology, number of constraints (NC) in 
the kinematic fit, and reconstruction efficiency, €, according 
to exclusive phase-space MC. 


Decay B [%] |2] Final state NC e [%] 
7 — yy 39.41+0.20 3y 1(n"r7) 5 26.58 
n — ntn m 22.92+0.28 3y 2(ntr7) 5 16.46 

Onr? 32.4840.23 Ty (rtr) 7 5.64 


7 — T 7 7 
Total 95.01+0.71 16.91 


xe obtained from four-momenta kinematic constraint 
fits are required to satisfy x20 < 40, xšce < 40 and 
x20 < 56 for n > yy, N > aha m) and n > 37°, 
respectively. These selection criteria effectively remove 
kaon and other charged track contamination, justifying 
the assumption that all charged tracks are pions. To 
select the Xcı candidates from the (3686) > yXcı tran- 
sition, we require the energy of the radiative photon to 
satisfy 0.155 < E, < 0.185 GeV. 


1. Background suppression 


The major background for all final states comes from 
(3686) > nJ/w, while in the n — yy case the back- 
ground from 7(3686) > yyJ/w decays is also significant. 
The background from (3686) — zzJ/4 is negligible, 
once a good 7 candidate is found. 


To suppress the (3686) > 7J/wW background for all 
three 7 decays, the system recoiling against the 7, with 
respect to the w(3686), must have its invariant mass sep- 
arated at least 20 MeV/c? from the J/w mass. 


Additional selection criteria are used in the n > yy 
channel to suppress 7° contamination and ~(3686) > 
yyJ/wv production. The former background is suppressed 
by rejecting events in which any two-photon combination 
satisfies 0.110 < m(yy) < 0.155 GeV/c?. The latter 
background is suppressed by vetoing events for which a 
two-photon combination not forming an 7 has a total 
energy between 0.52 GeV < E, < 0.60 GeV. This range 
of energies is associated with the doubly radiative decay 
(3686) > YXceJ; Xer > YJ /Y, for which the energy sum 
of two transitional photons is Hy, ~ 0.560 GeV. 


2. Background subtraction 


The background estimated from the inclusive MC after 
all selection criteria are applied is below 3% in each chan- 
nel. The background from 7) sidebands is subtracted, and 
Fig. shows the invariant mass distributions of 7 candi- 
dates with vertical dotted bars showing the 7 sideband 


regions. The sideband regions for the two-photon and 
three-pion modes are defined as 68 < |m(yy) — mn| < 
113 MeV/c? and 37 < |m(3m) — m,| < 62 MeV/c?, re- 
spectively, where m, is the nominal 7 mass B. In the 
case of 7 three-pion decays, the 7 signal region, defined 
by Eq. (i), is indicated by dash-dotted bars in Fig. 
Although the mass distribution of three neutral pions, 
Fig. [{c), is wider than the corresponding distribution 
from the charged channel, Fig. [i{b), we use the same 
selection criteria for both 7 decays, which keeps the ma- 
jority of good 7 — 37° candidates and results in similar 
background levels in the two channels. The effects of in- 
cluding more data from the tails of these distributions 
are taken into account in the systematic uncertainties. 
The invariant mass plot representing 7 — yy candidates, 
Fig. [I{a), is used only to select 7 sidebands for back- 
ground subtraction. Table [lists channel efficiencies and 
the effective efficiency for all channels. 


The nx*+7~ invariant mass distribution, when events 
from all 7 channels are combined, is shown in Fig. 2] In 
the signal region, indicated by vertical bars, there are 
33919 events, with the background of 497 events esti- 
mated from the 7 sidebands. The sideband background 
does not account for all the background, and after the 7- 
sideband background is subtracted, the remaining back- 
ground is estimated by fitting the invariant mass distri- 
bution. The fit is shown by the solid distribution, Fig. 
For the x< signal, a double-sided Crystal-Ball distribu- 
tion (dotted) is used, and for the background, a linear 
function along with a Gaussian corresponding to the x¢2 
contribution (dashed) are used. The signal purity esti- 
mated from the fit is P = (98.540.3)%, where the error is 
obtained from fluctuations in the background when using 
different fitting ranges and shapes of the background. 


B. Two-body structures in the xc — grt a 
decays 


The Dalitz plot for selected signal events is shown 
in Fig. Bla). Two-body structures reported in previous 
analyses of the xe > nwt decays, by BESII and 
CLEO [10, [27], the ao(980)z, a2(1320)m and f2(1270)n, 
are indicated by the long-dash-dotted, dashed and dash- 
dotted arrows pointing into the Dalitz space, respectively. 
One feature of this distribution is the excess of events in 
the upper left corner of the Dalitz plot (a), pointed to 
by the dotted arrows, which cannot be associated with 
known structures observed in previous analyses of this 
Xe1 decay. We hypothesize this is due to a2(1700) pro- 
duction. The expected Dalitz plot of a a2(1700)z signal is 
shown in Fig. Bb), obtained assuming that the a2(1700) 
is the only structure produced. The a2(1700) > nat 
and a2(1700) — 1 components cannot be easily iden- 
tified along the dotted arrows in the Dalitz plot, FigBla), 
but their crossing in the plot shown in Fig. Bb) visually 
matches the excess of events in the upper left corner of 
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FIG. 1. The invariant mass distribution of the 7 candidates, where dotted (red) lines indicate regions used for background 
subtraction, while dash-dotted bars (blue) show 7-signal boundaries for the three-pion 7 decay cases. There are no blue bars 
on plot (a) since the n — yy signal is selected using the yy kinematic constraint. 


the Dalitz plot of Fig. Bla). 


The distributions of the square of the invariant mass 
are shown in Fig. Blc) for nm and (d) for mta~. Struc- 
tures that correspond to ao(980), a2(1320) and f2(1270) 
production are evident, as well as a low-mass m7 peak, 
sometimes referred to as the o state. In each of these 
two distributions there is a visible threshold effect. In 
the zz distribution, there is a structure above the KK 
threshold, which is too broad to result from the fo(980) 
alone. In the ym distribution, the broadening of the 
ao(980) peak around 1.2 GeV*/c* could be associated 
with the n'r threshold. By examining various regions in 
the Dalitz space, we conclude that the cross-channel con- 
tamination, or reflections, are not associated with these 
threshold effects in the data. In order to eliminate back- 
ground as the source of these peculiar line shapes, back- 
ground studies are performed. Namely, we increased the 
background level by relaxing the kinematic constraint to 
XNc/NC < 10 and also suppressed more background by 
requiring y3¢/NC < 5. In addition we varied the limits 
on tagging 7 and xai candidates, as explained in Sec. 


It is possible that the zz line shape results from a 
destructive interference between the fo (980) and other 
components of the mz S-wave. It has been known for 
some time that the a6 (980) — nz line shape is affected 
by the proximity of the KK threshold to the ao(980) 
mass ba) If the a9(980) — n'r coupling appears to 
be important for describing the ag(980) — mz distribu- 
tion, this would be an example when a virtual channel 
is influencing the distribution of another decay channel, 
despite its threshold being far away from the resonance 
peak. We use an amplitude analysis (AA), described in 
the next section, to help in answering the above ques- 
tions, and to determine the nature and significance of 
the “crossing structure” discussed. 
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FIG. 2. Invariant mass of the x<i candidates, after the 


7 sideband background is subtracted. Vertical bars indicate 
the region used to select the x<i candidates. See the text for 
the fit discussion. 


II. AMPLITUDE ANALYSIS 


To study the substructures observed in the xa — 
mmm decays, we use the isobar model, in which it is 
assumed that the decay proceeds through a sequence of 
two-body decays, xa — Rhy; R — hihoə, where either 
an isospin-zero (R — Tm) or isospin-one (R — nr) res- 
onance is produced, with the total spin J, and relative 
orbital angular momentum L with respect to the bache- 
lor meson, hy. For resonances with J > 0, there are two 
possible values of L that satisfy the quantum number 
conservation for the 1++ — (JP°)0;7 transition. 


We use the extended maximum likelihood technique to 
find a set of amplitudes and their production coefficients 
that best describe the data. The method and complete 
description of amplitudes constructed using the helicity 
formalism are given in Ref. Bol, with two exceptions. 
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FIG. 3. Dalitz plots obtained from selected y-1 candidates from (a) data and (b) exclusive MC, assuming the a2(1700) is the 
only structure produced. The (c) 7m and (d) zT projections show various structures, which can also be identified by arrows 
in the Dalitz plot (a). Vertical dotted lines in plots (c) and (d) indicate the thresholds for producing the n'r or KK in the nr 


or 77 space, respectively. 


The first difference is that the events from the n- 
sidebands are subtracted in the likelihood function £, 
with equal weight given to the left-hand and right-hand 
sides, using a weighting factor w_= —0.5. The second 
difference with respect to Ref. io is that we deviate 
from the strict isobar model by allowing production am- 
plitudes to be complex. Isospin symmetry for nr® reso- 
nances is imposed. 


In the minimization process of the expression —2 ln £, 
the total amplitude intensity, I(x), constructed from the 
coherent sum of relevant amplitudes, is bound to the 
number of observed x-; candidates by using the integral 


Nas / €(x)I(x)dx, (2) 


where x represents the kinematic phase space, while €(x) 
is the acceptance function, with the value of one (zero) 
for accepted (rejected) exclusive MC events. The proper 
normalization of different 7 channels is ensured by using 
exclusive MC samples, generated with sample sizes pro- 
portional to the 7 branching fractions, listed in Table [I] 
If the complete generated exclusive MC set is used in the 
MC integration, then Eq. (2) provides the acceptance 
corrected number of xai events, adjusted by subtracted 
background contributions. In this case, (x) = 1 for all 
MC events. Fractional contributions, Fa, from specific 
amplitudes, Aa, are obtained by restricting the coherent 
sum in I(x) to a(x), so that 

f To(x)dx 


Fachan 


J I(x)dx (3) 
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The numerator represents acceptance-corrected yield of a 
given substructure, used to calculate relevant branching 
fractions, Ba. Errors are obtained from the covariance 
matrix using proper error propagation, so for a given 
substructure, the errors on Ba and Fa are not necessarily 
the same. 


The decay chain w(3686) + 7X1} Xe1 — NTT is de- 
scribed by amplitudes constructed to take into account 
the spin alignment of the initial state and the helicity of 
the radiated photon. Linear combinations of helicity am- 
plitudes can be used to construct amplitudes in the mul- 
tipole basis, matching the electric dipole (E1) and mag- 
netic quadrupole (M2) transitions. The (3686) > yx¢1 
decay is dominated by the £1 transition (CLEO) bo. 
and a small M2 contribution (~ 3%) can be treated as a 
systematic uncertainty. 


A. Mass dependent terms, T. (s) 


The dependence of amplitude Aa on the energy can 
be separated from its angular dependence, employing a 
general form p”q7T,(s), if the width of the xa is ne- 
glected. Here, p and q are decay momenta for decays 
Xa > Ryhy and R; > hiho in the rest frame of the Xe 
and a resonance RJ, respectively, while s = m7, is the 
squared invariant mass of the corresponding isobar prod- 
ucts (nr or nr). For most resonances, we use relativistic 
Breit-Wigner (BW) distributions, with spin-dependent 
Blatt-Weisskopf factors Bd. For the ao(980) and zz S- 
wave line shapes, we use different prescriptions explained 
below. 


To account for the nonresonant process Xel — NTT, 
we use an amplitude constructed as the sum of all pos- 
sible final state combinations of helicity amplitudes con- 
strained to have the same production strength, with no 
dependence on the invariant mass of the respective two- 
body combinations. 


1. Parametrization of ao(980) 


Instead of using the usual Flatté formula to de- 
scribe the a9(980) line shape, we use dispersion integrals, 
following the prescription given in Ref. . We consider 
three ag(980) decay channels, the nr, KK, and n'r, with 
corresponding coupling constants, geh, and use an appro- 
priate dispersion relation to avoid the problem of a false 
singularity present in the n'm mode (see the discus- 
sion at the end of this section). The ao(980) amplitude 
is constructed using the following denominator: 


Da(s) = må — s — 5 Ten(s), (4) 
ch 


where mo is the ap (980) mass and IL; (s) in the sum over 
channels is a complex function, with imaginary part 


Imll., (s) = JenPcn(8) Fen (s), (5) 


while real parts are given by principal value integrals, 
co / / 
Ralla eee L Di ie (6) 
T Je, (s-s) 

In the above expressions pep(s) is the available phase 
space for a given channel, obtained from the correspond- 
ing decay momentum qen(s): Pen(s) = 2den(s)/V/s. The 
integral in Eq. (6) is divergent when s > co, so the phase 
space is modified by a form factor F. (s) = e~2%n(s), 
where the parameter ñ is related to the root-mean- 
square (rms) size of an emitting source [20]. We use 
B = 2.0 [GeV /c?]~? corresponding to rms = 0.68 fm, and 
we verify that our results are not sensitive to the value of 
0. The integration in Eq. (6) starts from the threshold for 
a particular channel, sen, which conveniently solves the 
problem of the analytical continuation in special cases of 
final state configurations like the ao(980) — n'm, when 
the decay momentum below the threshold (s < my +m) 
becomes real again for s < My — Mr. Figure [4]shows the 
shapes of (a) ImII-n(s) and (b) Rellen (s), for the KK 
and mx channels, for arbitrary values of the coupling 
constants. In the final form, the real parts in the de- 
nominator of Eq. are adjusted by Rellen (mo) terms: 
Rellen(s) > Rellen(s) — Rellena (mo). 


2. nn S-wave model 


The zz S-wave parametrization follows the prescrip- 
tion given in Ref. , in which two independent pro- 
cesses for producing a mm pair are considered: direct 
(mm)s — (am)g, and production through kaon loops, 
(K K)s — (xm)g. Amplitudes corresponding to these 
scattering processes, labeled S,,(s) and Swgg(s), are 
based on di-pion phases and intensities obtained from 
scattering data , which cover the zz invariant mass 
region up to 2 GeV/c”. The S, (s) component is adapted 
in Ref. to account for differences in the m7 production 
through scattering and decay processes, using the denom- 
inator, D(s), extracted from scattering experiments. The 
S,(s) amplitude in this analysis takes the form: 


San (8) = coS*(s)+)) iZi er (S)S(8)+) | cz, (s)S°(s). 
i=l i=l (7) 


The common term in the above expression, S9(s) = 
1/D(s), is expanded using conformal transformations of 
the type 


VS + So — /Sth — 8 


z, (8) = ——— (8) 


Vs F so + /Sih — sS 


which is a complex function for s > sn. Equation (7) fea- 
tures two threshold functions, Zs„ (s), one corresponds to 
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FIG. 4. Line shapes of (a) ImII(s) and (b) Rell(s) for the KK and n'r production with arbitrary normalization. 


KK production with sgg = 4m%, while another with 
Sth = s’ could be used to examine other possible thresh- 
old effects in di-pion production. The c;, i = 1,2 are 
production coefficients to be determined. 


Figure[5] shows the (a) phase and (b) intensity of vari- 
ous components used in constructing the nm S-wave am- 
plitude based on two functions given by Eq. (8), with 
different thresholds: z¢g(s) and zs(s). The follow- 
ing convention is used: S$, (s) = zi gS (s), S (s) = 
z',$°(s). Components are arbitrarily scaled, and we set 
Vs’ ~ 1500 MeV/c’, similarly to the value used later in 
analysis. The parameter so = 1.5 (GeV/c?)? can be used 
to adjust the left-hand cut in the complex plane, and the 
same value is used in all components. 


IV. RESULTS 


We present results from the amplitude analysis of the 
full decay (3686) > YXc1; Yer > NTT, reconstructed 
in three major 7 decay modes. The optimal solution 
to describe the data is found by using amplitudes with 
fractional contributions larger than 0.5% and significance 
larger than 50. The significance for each amplitude a is 
determined from the change in likelihood with respect 
to the null hypothesis, AA = —2In£o/La. The null 
hypothesis for a given amplitude is found by excluding it 
from the base-line fit, and the corresponding amplitude 
significance is calculated taking into account the change 
in the number of degrees of freedom, which is two (four) 
for J = 0 (J > 0) amplitudes. 


The most dominant amplitude in this reaction is 
ao(980)7, as evident from the nr projection of the Dalitz 
plot, Fig. Bic). Other amplitudes used in our base- 
line fit include the SkRN, Sran, fo(1270)n, f4(2050)n, 
a2(1320)m and a2(1700)7, where masses and widths of 
resonances described by BW functions are taken from 
the PDG [9], while the a2(1700) and ao(980) parameters 


are free parameters to be determined by the fit in this 
work. The mass projections are shown in Fig. [6] and the 
corresponding fractional contributions and significances 
are listed in Table For amplitudes with spin J > 0 
both orbital momentum components are included. 


The following components form the S+, (s) amplitude: 
Sun (8) = o8. (5) + c1 Sar (8) + Cy Sem (s) +e STT (s). (9) 


As indicated earlier, the threshold used to construct the 
S*(s) term is spg = 4m%. The threshold for the S” (s) 
components (i = 1,2) is s’ = 2.23 [GeV/c?]?, which is 
close to the mass of the fo(1500), and it is responsible 
for the peaking of the S,,7 amplitude in this region, 
Fig.[6{b). In fact, the S”(s) components are used instead 
of the fo(1500)7 amplitude, which would be needed in the 
optimal solution if only threshold functions 2‘, = (s) were 
used in the expansion of the S,,(s)7 amplitude. With 
these additional terms, the contribution and significance 
of mm scalars, the fo(1370), fo(1500) and fo(1710), is 
negligible, for each. Although this particular set of am- 
plitudes respects the unitarity of the zz S-wave, we use 
the sum of BW to model other spins and final states, 
namely the f2(1270), f4(2050), a2(1320) and a2(1700). 
Our approach provides reasonable modeling of the zz 
line shape, and the sum of all mz S-wave components, 
Siz and S,,, is reported in Table [I] 


Besides the fo(1370), fo(1500), and fo(1710), other 
conventional resonances are probed, including the 
fo(1950), fo(1525), fo(2010), and a6 (1450), with param- 
eters fixed to PDG values a They do not pass the tests 
for significance and fractional contribution. The non res- 
onant xa — nwt 7 production is found to be negligible. 
The search for possible 1~* resonances in the mz final 
state will be presented below. 
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FIG. 5. The (a) phase and (b) intensity of the 7a S-wave components. Red (dashed) histograms represent the Sx amplitude, 


blue histograms (doted and dash-doted) are obtained using $ 


three-doted) represent S. = z SEa terms. 
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TABLE II. Fractional intensities F, and significances of amplitudes in the base-line fit, with the first and second errors being 
statistical and systematic, respectively. The third error for the branching fractions for the ye1 — Tr m decay and decays into 
significant conventional isobars is external (see text). For exotic mesons only statistical errors on their fractional contributions 
are provided. The upper limits for exotic meson candidates, which include both statistical and systematic uncertainties, are 
at the 90% confidence level. The coherent sum of all zz S-wave components, (zz )sn, is included in this report. Note, the 
branching fractions for amplitudes of the type Aan, involving isobars decaying into mtr”, are the products of ye1 > Aan and 


Aa — nr” rates. Branching fractions for isobars decaying into nr include charge conjugates. 
Decay Z [%| Significance [ø] B(xa > nnr ) [10-7] 
NTT T i3 4.67 + 0.03 + 0.23 + 0.16 
ao (980)* > 72.8 + 0.6 + 2.3 > 100 3.40 + 0.03 + 0.19 + 0.11 
a2(1320)" 7 3.8 + 0.2 + 0.3 32 0.18 + 0.01 + 0.02 + 0.01 
a2(1700)* 7 1.0:0.1 2 0.1 20 0.047 + 0.004 + 0.006 + 0.002 
Skk 2.5 = 0.2 = 0.8 22 0.119 + 0.007 + 0.015 + 0.004 
Sant 16.4 + 0.5 + 0.7 > 100 0.76 + 0.02 + 0.05 + 0.03 
(xtn7)sn 17.8 + 0.5 + 0.6 ú 0.83 + 0.02 + 0.05 + 0.03 
f2(1270)n 7.8 £0.34 1.1 > 100 0.36 + 0.01 + 0.06 + 0.01 
f4(2050)n 0.6 + 0.1 + 0.2 9.8 0.026 + 0.004 + 0.008 + 0.001 
Exotic candidates U.L. [90% C.L. 
74(1400) "7 0.58+0.20 3.5 < 0.046 
mı (1600) * 77 0.11£0.10 1.3 < 0.015 
m1(2015) * 77 0.060.03 2.6 < 0.008 


A. The a2(1700) signature 


All structures listed in Table [I] have been already re- 
ported in the decay y¥e1 — nnr, except the a2(1700)z. 
Its fractional contribution is around 1%, and the signifi- 
cance of each orbital momentum component is more than 
100. Detailed background studies are performed to en- 
sure that the background, remaining after 1-sideband 
subtraction, is not affecting the significance and frac- 
tional contribution of the a2(1700). Results of fitting 
the mass and width of the a2(1700), shown in Table [HI] 
are consistent with the values listed by the PDG l. To 


check how the a2(1700) parameters and fractional con- 
tributions are affected by the f2(1270) and a2(1320), we 
also fitted their masses and widths, which are provided 
in Table M] with statistical uncertainties only. The mass 
(width) of the f2(1270) is lower (higher) than its nominal 
value [2], maybe because of interference with underlying 
mm S-wave components or threshold effects, other than 
those for the KK or f9(1500) production. 


The systematic uncertainties for the a2(1700) mass and 
width are obtained by varying parameters of other am- 
plitudes within respective uncertainties listed in Ref. (ay, 
and taking into account variations listed in Table] The 
ao(980) errors are shown in Table Variations in the 
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a invariant mass from data, compared with our base-line fit (solid curve) and 


corresponding amplitudes (various dashed and dotted lines). All features of the data, including structures discussed in Sec. 


are reproduced rather well. 


shape of the 7a-S wave amplitude are taken into account 
by changing terms in the expansion, Eq. (9). 


We also test the significance of the a2(1700) includ- 
ing alternative states with the same mass and width, but 
different spins: J = 0,1,4. In all cases, the significance 
of the a2(1700) in the presence of an alternative state 
exceeds 17a. The statistical significance of the a2(1700) 
signal alone is 200. This result confirms our hypothesis 
based on a visual inspection of the Dalitz plot, Fig. Bla), 
that the excess of events in the upper left corner of the 
Dalitz space results from the a2(1700) production, and 
it is associated with the crossing of the a2(1700)ta— 
and a2(1700)~ z" components. Further, Fig. [shows the 
nm mass distribution in the region around the expected 
a2(1700) peak, where data points are compared with a 
fit when the a2(1700)z amplitude is excluded. 


B. ao(980) parameters 


When determining the ao (980) parameters we use the 
ratios Rai = Ge g/ Gans and R31 = gy gb. + The result- 
ing values are listed in Table[[V] where systematic uncer- 
tainties are obtained by fitting the a9(980) parameters 
under different conditions. The level of background is 
varied by changing selection criteria described in Sec. M 
and by changing the amount of background subtracted 
from the 7 sidebands. Effects of the line shapes of the 
a2(1320), a2(1700), f2(1270) and f4(2050) resonances are 
taken into account by varying their masses and widths 
within the respective uncertainties Bl. and using values 
from Table [I] The effect of the zz S-wave shape is ex- 
amined in similar way as for the a2(1700). The presence 
of alternative conventional and exotic resonances is also 
taken into account. Our result is not sensitive to the 
value of the parameter 8 in Eqs. (5) and (6), within the 
range of values: 8 = (2.0 + 1.0) [GeV/c?]?. 


For comparison we list two previous results, one from 
a similar experiment, CLEO-c, and the other obtained 


using Crystal Barrel data. There is a general agreement 
between different analyses for the ag(980) mass and Rar. 
The ratio R3; was fixed in Ref. to the theoretical 
value provided by Eq.(11), while it was consistent with 
zero in the CLEO-c analysis, possibly because of smaller 
statistics. It is not easy to comment on the difference in 
values for the nz coupling, which could be affected by 
different normalizations used by different analyses. 


This analysis provides the first nonzero measurement 
of the coupling constant gwr. To test the sensitivity of 
the a6 (980) — nz line shape to the decay ag (980) > mz, 
we repeat the analysis with gyr = 0, and let the values 
of the other parameters free. The results of this fit are 
also given in Table The likelihood change when the 
n'm channel is ignored shows that the significance of a 
nonzero Gy’x Measurement is 8.97. The same result is 
obtained when the analysis is performed in the presence 
of the a6 (1450). The values of the two ratios based on 
the SU(3) expectation are 


g% k /g2, = 1/(2cos? $) = 0.886 + 0.034, 


(10) 


gya | Ion = tan? ó = 0.772 + 0.068, (11) 
which depend on the choice of the 7—1 mixing angle; ¢ = 
(41.341.2)° in this case [20]. Our result is consistent with 
Eq. within 1.50, based on the quoted uncertainties. 


C. Search for nz P-wave states 


We examine possible exotic meson production in the 
nm invariant mass region from 1.4 to 2.0 GeV/c?. Ta- 
ble [LI] lists fractional contributions and significances of 
three JP? = 1~-+ candidates, added one at the time to 
our nominal fit. Two possible orbital-momentum config- 
urations for an exotic amplitude are the S-wave and D- 
wave, and the significance of each is tested individually. 
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FIG. 7. 


The nz invariant mass projection from data in the region (1.3, 2.4) GeV/ c, compared with the fit without the 


a2(1700)7 amplitude (solid curve). Other amplitudes are plotted (various dashed and dotted lines) for comparison, while the 


peak that is associated with the a2(1700) is evident. 


TABLE III. The mass and width of the a2(1700), with statistical and systematic uncertainties. Only statistical uncertainties 
from the f2(1270) and a2(1320) fits are listed. Comparison with the PDG [A values is provided, with all units in GeV/c’. 


BESIII PDG [2] 
Resonance M r M I 
a2(1700) 1.726+0.012+0.025 0.190+0.018+0.030 1.732+0.016 0.194+0.040 
fə (1270) 1.258+0.003 0.206+0.008 1.275+0.001 0.185+0.003 
a2(1320) 1.317+0.002 0.090+0.005 1.3818+0.001 0.107+0.005 


We find that the significance of the S-wave is marginal, 
less than 2o for every 7, and the reported significances 
in Table M]result from using the S and D waves together 
in the fit. The most significant of the three possible ex- 
otic states is the ri (1400), with a significance of 3.50 and 
fractional contribution less than 0.6%. This represents a 
weak evidence for the existence of the 71(1400) because 
in alternative amplitude configurations, when parame- 
ters of other amplitudes are varied, the significance of 
this state becomes < 30. In the nominal amplitude con- 
figuration, the significance of each zi (1400) component is 
less than 30, and when taken together, the contribution 
of the S-wave is much smaller than the D-wave contri- 
bution, pointing that the evidence for the 71(1400) is 
circumstantial. 


Masses and widths of the three exotic candidates are 
not very well constrained by previous analyses, and we 
vary the respective parameters within listed limits 
Our conclusion is that there is no significant evidence for 
an exotic nm structure in the ye. > nrn decays, and 
we determine upper limits at the 90% confidence level for 
the production of each mı candidate. 


D. Branching fractions 


The branching fraction for the yaa — nwt a7 decay is 
given by 


P x N. = 
B(x > natn) = it 


= —P*N eon 49 
Nop(3686) By (3686) 37x01 Pn€ 


where the branching fractions Basse) yxa. and By are 
from Ref. PË the latter is listed in Table [I] The num- 
ber of (3686), Nusese); is provided in Sec [Tl] 
The signal purity, P, given in Sec. takes into ac- 
count that the number of xai obtained from the ampli- 
tude analysis includes the background not accounted for 
by the sideband subtraction. Using Eq. (2) we obtain 
Ny: = 192658 + 1075, where the error is from the co- 
variance matrix. The efficiency in Eq. is € = 1, by 
construction. 


Table [LI] lists the branching fraction for the xe — 
nnr, and branching fractions for subsequent reso- 
nance production in respective isospin states, mz or 
m+n, where the first and second errors are statistical 
and systematic, respectively. The branching fraction for 
a given substructure is effectively a product: 


Ba = Fa x B(Xe1 =? mm), (13) 


obtained using generated exclusive MC in accordance 
with Eq. (3). The third error is external, associated with 


201711.02325v1 


chinaXiv 


13 


TABLE IV. Parameters of the ao (980) determined from the fit using the dispersion relation of Eqs. (416), compared to results 
from previous analyses. Bold values indicate quantities that are fixed in the fit. 


Data mo [GeV/c'] Inr [GeV/c] 9k Gne gaz Inz 
CLEO-c [10] 0.998 + 0.016 0.36 + 0.04 0.872 + 0.148 0.00 + 0.17 
C.Barrel [20] 0.987 + 0.004 0.164 + 0.011 1.05 + 0.09 0.772 
BESIII 0.996+0.002+0.007 0.368+0.003+0.013 0.931+0.028+0.090 0.489+0.046+0.105 
BESIII (Rš, = 0) 0.990+0.001 0.341+0.004 0.892+0.022 0.0 


uncertainties in the branching fractions for the radiative 
transition (3686) > yxe1 and 7 decays. We also show 
the total rtn S-wave contribution, obtained from the 
coherent sum of the Sgg and S;, components. Statis- 
tical errors, as well as systematic ones, for a given frac- 
tional contribution and branching fraction differ, because 
common systematic uncertainties for all amplitudes can- 
cel when fractions are calculated, which will be discussed 
below. 


The upper limits for the production of the 7 (1400)z, 
71 (1600), and zi (2015)z are shown in Table] The lim- 
its are determined by including the corresponding ampli- 
tude in the nominal fit, one at a time. The analysis is 
repeated by changing other amplitude line shapes, and 
the background level, in a similar fashion used for deter- 
mining systematic uncertainties of nominal amplitudes 
(see Sec. [V). Masses and widths of exotic candidates are 
also varied within limits provided by the PDG B. The 
largest positive deviation of the exotic candidate yield 
with respect to the corresponding yield from the mod- 
ified nominal fit is effectively treated as the systematic 
error, summed in quadrature with the statistical error 
on a given exotic state yield. The resulting uncertainty 
is used to determine the 90% confidence level deviation, 
and added to the ”nominal” yield of an exotic candidate 
to obtain the corresponding upper limit for the branching 
fraction B(X > TITT). 


The branching fractions for the substructures in x= > 
mmm decays reported by the PDG |2] are compared in 
Table[V] with the values measured in this work, and with 
the previous most precise measurement (CLEO-c) id. 
The measurement for the f2(1270) production is adjusted 
to account for the measured relative f2(1270) > ntr- 
width. There is a rather large discrepancy between the 
values for the two most dominant substructures listed 
by the PDG and the two most recent measurements. 
There is very good agreement between the last two mea- 
surements, suggesting that the PDG values on two-body 
structures observed in Xe1 — nnt rT need to be updated. 


V. SYSTEMATIC UNCERTAINTIES 


Tables[VI]summarizes various contributions to the sys- 
tematic uncertainties in determining the Xe — Tr! 


TABLE V. Comparison between recent measurements of the 
branching fractions B(ye1 — natn), and with the PDG 
values. 


B(Xa — mr 7) x [10 2] 
Decay BESIII CLEO-c [10] PDG [2] 
Tr T 4.67 + 0.28 4.97 + 0.31 4.9 + 0.5 
ao(980)" x 3.40 + 0.23 3.29 + 0.22 1.8 0.6 
f2(1270)n 0.64: 0.11 066+011 2.7+0.8 


branching fraction, and Table shows the systemat- 
ics on the fractional contributions of amplitudes in the 
nominal fit. Systematic uncertainties in determining the 
Xa — mrt branching fraction stem from uncertain- 
ties in charged track and shower reconstruction efficien- 
cies, the contribution of the M2 multipole transition, 
amplitude modeling, the background contribution, and 
the uncertainty in the number of (3686) produced at 
BESHI (23) . External sources of uncertainty include 
the branching fraction B(w(3686) + 7x1) and the frac- 
tion of 7 decays, B(n) in Eq. (12). The external error 
affects only branching fractions, not fractional contribu- 
tions, and it is reported as a separate uncertainty. 


Systematic uncertainties associated with the tracking 
efficiency and shower reconstruction are 1% per track and 
1% per photon. Because of different final states used 
in this analysis, tracking and photon uncertainties are 
weighted according to the product of branching fractions 
and efficiencies of the different 7 channels, as listed in Ta- 
ble [I] The resulting systematic uncertainties for charged 
tracks and photons are 2.47% and 3.92%, respectively. 


The electromagnetic transition w(3686) — Xa is 
dominated by the £1 multipole amplitude with a small 
fraction of the M2 transition g). The nominal fit takes 
only the E1 multipole amplitude. Adding a small con- 
tribution of the M2 helicity amplitude, of 2.9%, we find 
a difference in the branching fraction of 0.62%. This is 
taken as a systematic uncertainty. 


When considering the effects of modeling line shapes 
of different amplitudes, we repeat the analysis changing 
the mass and width of resonances, a2(1320), f2(1270), 
and f4(2050), within respective uncertainties, and change 
the ao(980) and a2(1700) parameters within the lim- 
its of their statistical uncertainties, given in Tables 
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and We also change BW line shapes by replacing 
spin-dependent widths with fixed widths, and take into 
account the x<i width and centrifugal barrier as another 
systematic error. The largest effect from all these sources 
is taken as a systematic uncertainty for the branching 
fractions and fractional contributions. 


The effect of background is estimated by varying 
the kinematic-constraint requirement, changing limits on 
tagging 7 and xc<i candidates, changing the level of sup- 
pression of the J/7 and 7°? productions, and the level 
of background subtraction. As a general rule, selection 
criteria were changed to allow for ~ lo additional back- 
ground events, based on the numbers from the inclusive 
MC. We use xwe/NC < 9 in all three modes when 
varying the kinematic constraint. Based on these vari- 
ations, we conclude that the systematic uncertainty as- 
sociated with the assumption that all charged tracks are 
pions is negligible. To select x<i candidates, we use pho- 
ton energy ranges of (0.152-0.187) GeV, in the n > yy 
channel, and (0.150-0.190) GeV, in two 7 — 37 chan- 
nels. The mass window for the 7 selection is changed to 
(0.530-0.565) GeV /c?. The 7° suppression window is re- 
duced to (0.120-0.150) GeV/c? and the J/~ suppression 
is reduced by vetoing two-photon energy within (0.525- 
0.595) GeV. We also determine the branching fractions 
without background subtraction from -sidebands, and 
the largest effect is listed in Tables [VI] and [VI] 


Some uncertainties that are common for all ampli- 
tudes, like tracking, shower reconstruction, and N,,(3686) 
errors, cancel out in the fractional contributions. How- 
ever, they are taken into account when branching frac- 
tions are determined. 


TABLE VI. Systematic uncertainties in determining the 
branching fraction B(xe1 — natn). The systematic un- 
certainty per track is 1.0%, and for photons it is 1.0% per 
shower. 


Contribution Relative uncertainty (%) 
MDC tracking 2.5 
Photon detection 3.9 
M2/E1 0.6 
Background 1.6 
Amplitude modeling 0.1 
N. ~ (3686) 0.7 
Total 5.0 
External 3.4 


VI. SUMMARY 


We analyze the world’s largest Xe1 — mr! 7 sample, 
selected with very high purity, and find a very promi- 
nent ao(980) peak in the n7* invariant mass distribution. 
An amplitude analysis of the 7(3686) > yxe1; Xea > 
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TABLE VII. Systematic uncertainties in fractional contribu- 
tions, in percent, for the base-line amplitudes used to model 
the xa 3 nwt a decays. 


Source M2/E1 Background T,(s) Total 
ao (980) T 0.2 05 3.1 32 
a2(1320)r 0.5 5.6 5.6 7.9 
a2(1700)r 1.4 3.8 12 13 
SkkN 3.7 2:2 11 11.5 
Sppn 1.1 li 4.3 4.6 
TTSN 1.5 1.1 3.0 3.6 
f2(1270)n 0.5 2.3 14 15 
Ja (2050) 5.6 25 18 32 


natm decay is performed, and the parameters of the 
ao (980) are determined using a dispersion relation. The 
ao(980) line shape in its mz final state appears to be sensi- 
tive to the details of the a9(980) — mz production, and 
for the first time, a significant nonzero coupling of the 
ao(980) to the n'm mode is measured with a statistical 
significance of 8.90. 


We also report a2(1700)z production in the xa > 
natn decays for the first time, with the mass and width 
in agreement with world average values, and this analy- 
sis provides both qualitative and quantitative evidence 
for the existence of the a2(1700). First, the signature of 
the a2(1700) in the Dalitz space is consistent with the ob- 
served Dalitz plot distribution. Second, the a2(1700) sig- 
nificance from the amplitude analysis is larger than 17<, 
compared to alternative spin assignments, even though 
the fractional yield of the a2(1700)m is only 1%. This 
may help in listing the a2(1700) as an established reso- 
nance by the the PDG Jl. 


We examine the production of exotic mesons that 
might be expected in the xa — nam decays: the 
71(1400), 71(1600) and 7(2015). There is only weak 
evidence for the 7(1400) while other exotic candidates 
are not significant, and we determine the upper limits on 
the respective branching fractions. 
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